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, We study the effects of pressure anisotropy and heat dissipation in a spherically symmet 

C "3 ■ ric radiating star undergoing gravitational collapse. An exact solution of the Einstein 

field equations is presented in which the model has a Friedmann-like limit when the 
heat flux vanishes. The behaviour of the temperature profile of the evolving star is in- 
vestigated within the framework of causal thermodynamics. In particular, we show that 
$H ■ there are significant differences between the relaxation time for the heat flux and the 

relaxation time for the shear stress. 
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1. Introduction 



It has been shown that the role of anisotropy can alter the evolution and subse- 
quently the physical properties of stellar objects. As an example, investigations have 
shown that the maximal surface redshift for anisotropic stars may differ drastically 
from isotropic stars. The origin of anisotropy within the stellar core has received 
widespread attention amongst astrophysicists. A review of the origins and effects 
of local anisotropy in stellar objects was carried out by Herrera and Santos 1 , and 
Chan et aP and more recently by Herrera et aP. The physical processes that are 
responsible for deviations from isotropy can be investigated in the high and low 
density regimes. Hartle et aft have shown that pion condensation at nuclear densi- 
ties (0.2/ -3 < p < 2.0/~ 3 , If = 10~ 13 cm) can drastically reduce the pressure and 
hence impact on the evolution of the collapsing star. At higher densities the short 
range repulsion effects dominate which damp out the pionic effects giving rise to 
significantly different values for the pressure. As pointed out by Martinez 5 viscos- 
ity effects due to neutrino trapping at nuclear densities can alter the gravitational 
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collapse of a radiating, viscous star. Anisotropic velocity distributions and rotation 
can induce local anisotropy in low density systems. 

The boundary of a radiating star divides spacetime into two distinct regions, 
the interior and the exterior region. In order to fully describe the evolution of such a 
system, one needs to satisfy the junctions conditions for the smooth matching of the 
interior and exterior spacetimes. In the case of a shear-free, spherically symmetric 
star undergoing dissipative gravitational collapse these conditions were first derived 
by Santos 6 . Subsequently, many models of shear- free radiative collapse were devel- 
oped and the physical viability of these models were studied in great detail 7 ' 8 ' 9,10,11 . 
These junction conditions were later generalised to include pressure anisotropy 12 
and the electromagnetic field 13 . 

There are a few exact solutions to the Einstein field equations for a bounded 
shearing matter configuration. There have been numerous attempts to pro- 
duce models of radiative gravitational collapse which incorporate the effects of 
shear 14 ' 15 ' 16 . Most treatments to date are based on numerical results as the re- 
sulting temporal evolution equation derived from the junction conditions is highly 
nonlinear 15 . It is in this spirit that we seek an exact solution of the Einstein field 
equations which represents a spherically symmetric radiating star undergoing dis- 
sipative gravitational collapse with nonzero shear. 

In this treatment, we consider the relaxation effects due to the heat flux and 
shear separately. We show that earlier assumptions of the relaxation time being 
proportional to the corresponding collision time only hold for a limited regime of 
the evolution of the star. These results agree with earlier suggestions by Anile et 
al . We are also in a position to integrate the causal heat transport equation and 
obtain the corresponding causal temperature profile in the interior of the star. 

This paper is organised as follows. In section 2, we provide the Einstein field 
equations for the most general, nonrotating, spherically symmetric line element. 
The energy momentum tensor for the interior spacetime is that of an imperfect fluid 
with heat conduction and pressure anisotropy. In section 3, the junction conditions 
required for the smooth matching of the interior and exterior spacetimes across a 
time-like hypersurface are presented. The transport equations for the dissipative 
fluxes within the framework of causal thermodynamics are presented in section 4. 
The general solution for the special case of constant mean collision time in both 
the causal and noncausal theories are presented. The junction conditions are used 
to generate a reasonable model of radiating anisotropic collapse in section 5. When 
the heat flux vanishes, our model tends towards a Friedmann-likc limit, similar to 
the result obtained for the shear-free case studied by Kolassis et aF . A discussion 
of some relevant physical parameters is presented in section 6. 
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2. Spherically symmetric Spacetimes 

The interior spacetime of a non-rotating spherically symmetric collapsing star is 
most generally described by the line element (in coordinates (x a ) — [t,r,6,4>)) 

ds 2 = -A 2 dt 2 + B 2 dr 2 + Y 2 {d9 2 + sin 2 9d0 2 ) , (1) 

where A, B and Y are functions of the coordinates t and r. The fluid four-velocity 
u is comoving and is given by 

«° = ^<5o- (2) 
The kinematical quantities for the line element (1) are given by 

°4(§+4). 

where u a is the four-acceleration vector and O is the expansion scalar. The nonzero 
components of the shearing tensor are given by 2 

'/ \ 

(5) 




(6) 

033 = C22 sin 2 6. (7) 
The magnitude of the shear scalar is given by 

'-^-'-£s(|-f)- (8 » 

For a relativistic fluid, the kinematical quantities are important for studying the 
evolution of the system. The interior matter distribution is described by the follow- 
ing energy-momentum tensor 

Tab = (P+ p)u a Ub + P9ab + TTah + q a U b + q b U a , (9) 

where p is the isotropic pressure, p is the density of the fluid, Tr a b is the stress tensor 
and q a is the heat flux vector. The stress tensor takes the form 

TTab = (PR - PT)(n a n b - ^h ab ) , (10) 

where pn is the radial pressure, pr is the tangential pressure, and n is a unit radial 
vector given by 

«° = ^. (11) 
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The isotropic pressure p is related to the radial pressure and the tangential pressure 
via 

p=^[ PR + 2 PT }. (12) 
The coupled Einstein field equations for the interior matter distribution become 

2 BY 1 1 Y 2 
P ~ A 2 BY + Y 2 + A^Y 2 

1 / Y" Y' 2 B'Y'\ 

-w( 2 - + W- 2 bT) < 13 > 

1 ( Y Y 2 AY\ 

Pr = T 2 \- 2 y-- 2+2 ay) 

1 / Y' 2 A' Y'\ 1 

+bi{y> +2 ay)-y- 2 (14) 

1 / B AB BY AY y\ 
Pt ~ \ B~AB + by ~ay + y) 




'/ -JJjl I -- --)• ^ ! «' 

In the above we have defined 

q = q\ (17) 

where q a = (0, q 1 ,^), 0). The field equations (13)-(16) describe the gravitational 
interaction of a shearing matter distribution with heat flux and anisotropic pressure. 

3. Spherical collapse with heat flow 

The problem of gravitational collapse within the context of general relativity was 
first investigated by Oppenheimer and Snyder 18 in which the interior matter dis- 
tribution was taken to be dust with the exterior spacetime being Schwarzschild. 
With the discovery of the Vaidya solution 19 , it became possible to study radia- 
tive gravitational collapse where the collapsing core radiated energy to the exterior 
spacetime 6 ' 8, 16,20 . The Vaidya solution which describes the exterior spacetime of a 
radiating star is given by 

ds 2 = - (l - dv z _ 2dvdr + r 2 (d9 2 + sin 2 9d<f> 2 ) (18) 
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in the coordinates (x a ) = (v,r,6,(j)). The quantity m(v) represents the Newtonian 
mass of the gravitating body as measured by an observer at infinity. The metric 
(18) is the unique spherically symmetric solution of the Einstein field equations for 
radiation in the form of a null fluid. The Einstein tensor for the line element (18) 
is given by 

G °»=-^S>°- (i9) 

The energy-momentum tensor for null radiation assumes the form 

Tab = ew a w b , (20) 

where the null four-vector w is given by w a = (1, 0, 0, 0). Thus from (19) and (20) 
we have that 

2 dm . . 

^-7^ (21) 

for the energy density of the null radiation. Since the star is radiating energy to the 

dm 

exterior spacctimc we must have — — < 0. 

dv 

In deriving the junction conditions we employ the Darmois conditions as these 
were shown to be the most convenient and reliable 21 . Here we merely state the 
results of the matching of the line elements (1) and (18) since these conditions have 
been extensively investigated in the past. For a more comprehensive treatment of 
the junction conditions the reader is referred to the works of Glass 22 , Chan 12,15 and 
Govender et al 11 . 

We require that the metrics (1) and (18) match smoothly across the boundary 
E. This generates the first junction condition 

(ds 2 _h = (<&+)e =ds%. (22) 

(We use the notation ( )s to represent the value of ( ) on £.) The second junction 
condition is obtained by requiring continuity of the extrinsic curvature across the 
boundary. This gives 



*S = *«. (23) 



where 



+ + d 2 Xl , „ dXldXl , N 

X'^-^We^-^WW- (24) 

and n„ (X±) are the components of the vector normal to S. We find that the 
necessary and sufficient conditions on the spacetimes for the first junction condition 
to be valid are that 

i 

A{r s ,t)dt= (l-^ + 2^) a dv (25) 



Y{rv,t) = rv{v). (26) 
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By equating the appropriate extrinsic curvature components for the interior and 
exterior spacetimes we generate the second set of junction conditions which arc 
given by 



m(v) = 



1 



A 2 



Y 2 

1 r 



(27) 



(Pr)e = (g-B) E , 



(28) 



where we may interpret m(v) as representing the total gravitational mass within 
the surface S. The expression (27) corresponds to the mass function of Cahill and 
McVittie representing spheres of radius r inside S. 

The important result (pr)s — (qB)s, relating the radial pressure pr to the heat 
flow q, was first established by Santos 6 for shear-free spacetimes. The first attempt 
to generalise the above junction conditions to include shear for neutral matter was 
carried out by Glass 22 . 

The equations (25)-(28) are the most general matching conditions for the spheri- 
cally symmetric spacetimes A4 + and Ai~ . Relation (28) implies that the radial pres- 
sure pr is proportional to the magnitude of the heat flow q which is non-vanishing 
in general. 

The total luminosity for an observer at rest at infinity is given by 



,W " 1 dv)~ 2 



Y 



+ Y 



(29) 



dm 



where < since is positive. An observer with four-velocity v a = (v, r, 0, 0) 

dv 

located on E has proper time r\ related to the time t by d-q — Adt. The radiation 
energy density that this observer measures on S is 



v 2 dm 



dv 



and the luminosity observed on S can be written as 

L E = 47rr 2 e E . 

The boundary rcdshift z E of the radiation emitted by the star is given by 

-l 



1 + zs 



dv 
d-q 



~B 



Y 



(30) 



(31) 



(32) 



which can be used to determine the time of formation of the horizon. The above 
expressions allow us to write 



is 



(33) 
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which relates the luminosities L^, to L M . The redshift for an observer at infinity 
diverges at the time of formation of the horizon which is determined from 

Y' 



B 



+ Y = 0. 



(34) 



4. Thermodynamics 

Previous treatments of heat flow in relativistic stellar models have shown that 
relaxational effects play a significant role in the evolution of the temperature and 
luminosity profiles during late stages of collapse 5 ' 10 ' 11 ' 24 . Since we are interested in 
the effects of shear within the stellar core, we employ the causal transport equations 
for the heat flux and the shear stress. Assuming there is no viscous/heat coupling, 
we have the following relevant transport equations for the dissipative fluxes 25 



Th a b q b + q a = -n(D a T + Tu a ) - 



i r 2 



T\h a c h h -k cd + n ab = -2t]a ab 



\kT 2 



T^ab 



q„ 



(35) 



(36) 



where r and t\ are the relaxation times of thermal and shear viscous signals re- 
spectively. The truncated transport equations, together with the no-coupling as- 
sumption are differential equations of Maxwell-Catteneo form 

Th a b q b + q a = —n(D a T + Tu a ) (37) 



T\h a c h b d 'k c d + TT ab — -2r\o ab . 



(38) 



Setting r = t\ = in the above, we regain the so-called Eckart transport equations 
which predict infinite propagation velocities for the dissipative fluxes. 

In this paper we are primarily interested in heat transport in relativistic as- 
trophysics and hence (37) plays a significant role in determining the evolution of 
the causal temperature profile of our models. For the line clement (1) the causal 
transport equation (37) becomes 



T(qB) tt + A(qB) = -k- 



B 



which governs the behaviour of the temperature. Setting r 
the familiar Fourier heat transport equation 

{AT),r 



A(qB) 



B 



(39) 

in (39) we obtain 
(40) 



which predicts reasonable temperatures when the fluid is close to quasi-stationary 
equilibrium. 

For a physically reasonable model, we use the thermodynamic coefficients for 
radiative transfer outlined in Martinez 5 . We consider the situation where energy is 
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carried away from the stellar core by massless particles that are thermally generated 
with energies of the order of kT. The thermal conductivity takes the form 

« = 7T 3 T C , (41) 

where 7 (> 0) is a constant and r c is the mean collision time between the massless 
and massive particles. Based on this treatment we assume the power-law behaviour 

r c = (2) T-» , (42) 

where a (> 0) and w (> 0) are constants. With lo = | we regain the case of 
thermally generated neutrinos in neutron stars. The mean collision time decreases 
with growing temperature, as expected, except for the special case to = 0, when 
it is constant. This special case can only give a reasonable model for a limited 
range of temperature. Following Martinez 5 , we assume that the velocity of thermal 
dissipative signals is comparable to the adiabatic sound speed which is satisfied if 
the relaxation time is proportional to the collision time: 

'/V 



Q 



(43) 



where r (> 0)is a constant. We can think of r as the 'causality' index, measuring 
the strength of relaxational effects, with r = giving the noncausal case. 
Using the above definitions for r and n, (39) takes the form 

0{qB) tt T-" + A(qB) = -a ^-^ . (44) 

When f3 = 0, we can find all noncausal solutions of (44), viz. 

(AT)*'" = ^— l - [ A 4 ~"qB 2 dr + F(t) lo ± 4 
a J 

In (AT) = -- [ qB 2 dr + F(t) uj = 4, (45) 
a J 

where F(t) is an arbitrary function of integration. This is fixed by the expression 
for the temperature of the star at its surface S. 

In the case of constant mean collision time, ie. lo — 0, the causal transport 
equation (44) is simply integrated to yield 

4 " 



(AT) 



J A 3 B(qB), t dr + J A 4 qB 2 dr 



+ F(t) (46) 



while one solution valid for a less limited range of temperature can be found for 
ui = 4, which corresponds to nonconstant collision time 26 : 

(AT) 4 = -M exp (-/^dr) J A 3 B(qB) it exp (/ ^dr) dr 

+ F(t) exp (- j ^^dr") . (47) 
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In order to investigate the relaxational effects due to shear we utilise (38) as a 
definition for the relaxation time for the shear stress. For the metric (1) the shear 
transport equation (38) reduces to 



-P 



T\ = 



P+±r oT*' 



(48) 



where we have used the coefficient of shear viscosity for a radiative fluid 25 

P = 3 (pt — Pr) and r is the radiation constant for photons. We have further 
assumed that t\ = &\t c . 



(49) 



5. Radiating anisotropic collapse 

The junction condition (pr)s = (<7-B)s yields 



2YY + Y 



>2 



Y" 2 
B 2 " + B 



D 



YY' - 2—YY' + 1 
B z 



0, 



where we have set A(r, t) = 1. A particular solution of (50) is given by 



it 1 / 3 



B 

which yields the line clement 
ds 2 = -dt 2 + 



1 + Ci(r)e 
, 1 - ci(r)e S1 ^ 



1 + ci(r)e~ 



f 2/3 



dr 2 + r W 



,1 - ci(r)e r ) 

With this form of the line element, the field equations (13)-(16) become 
/ / 



P = 



3i 2 



1 + 



r 2 U 



18 1 



\ 



l-e^ci(r) (l + e^c^r) 



3(rt§+9i) - (r 2 + 3 r t'i + 9 1 



1 + e r ci(r) 



1 + e r ci(r) 



3e— t3 -1 + e— ci(r) ci(r) 1 



r I 1 + e 3 * 3 ci (r) 



(50) 

(51) 
(52) 

(53) 



(54) 
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PR = 



PT 



— 4e r ci(r) 



t$ I r + e 3t ^ rc\{r) 



3r 3 *s 



-3 r 1 3 



+ 



l + e^-ci(r) -1 + e— ci(r) 



9*3 



+ 



3 (3*3 ci(r) - r 2 c'^r)) 6*3 (-3*3 ci(r) + r 2 ci(r)) 



ci(r) ( 1 + e 3 *' ci(r) 



ci(r) ( 1 + e 3 * 3 ci(r) 



2 r 2 ci (r) - 9 * S ci (r) + 3 r 2 * 3 ci (r) 
I 

/ \ . 3*3. , ,2 

ci(r) + e - ci(r) 



+ 



3 = 



3t? / \ I , 3*5 , , 

4e »• ci(r) —1 + e r Ci(r) 



r 2 * 2 l + e^ciW 



In the absence of heat flux (ci = 0) our model yields 

Y = rt 2 '* 



(55) 



(56) 



(57) 



(58) 



B = * 2 / 3 (59) 

9 = ± (60) 

Pr=Pt = 0. (61) 

(Note that the solution presented above is not simply obtained by substituting 
ci = into (54)-(57). Rather, one must go back to the original metric and then 
derive the field equations ab initio.) The above solution represents a dust sphere 
and the metric is described by the Einstcin-de Sitter solution. We note that when 
q = the pressure must vanish which allows for the matter to have free-fall motion. 
For q ^ 0, the pressure is non-vanishing so that it compensates for the outgoing 
heat flux thus allowing for free-fall motion. With this in mind we expect that the 
luminosity radius of the star in both the radiative and non-radiative cases have 
the same temporal dependence. Calculating the luminosity radius for our radiating 
model, we obtain 



(62) 
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Fig. 1. Causal temperature (dashed line), noncausal temperature (solid line) versus r. 



which is independent of ci . Hence the case ci = reduces to the model investigated 
by Oppcnhcimcr and Snyder 18 . 



6. Physical considerations 

In order to check the physical viability of our model, we investigate the evolution 
of the temperature profile within the framework of extended irreversible thermody- 
namics. 

Utilising (46) and (53) we obtain 



16/3ci 

(47T(5r 2 ) E ' 3crf 5 / 3 



+ 



a 3t L 



6(l + cie 3tl/3 /6) r(l + cie 3 * 1/3 A-) 



16 



log 



-1 + cie 



3t 1/3 /6" 



1 + cie 



3t' 



-1 + Cl e 3tl/3 /- ) \l + c 1 e 3tl/3 / b 



+ 



16 

Sat 



tanh-^cie 24 /h ) - tanh _1 (cie 



2t'' A /r^ 



(63) 



where is given by (29). We note that the causal and the noncausal ((3 = 0) 
temperatures coincide at the boundary (r = b): 

T(t,b) = f(t,b) . (64) 

However, Figure 1 shows that at all interior points, the causal and non-causal tem- 
peratures differ. In particular, we observe that the causal temperature is greater 
than the non-causal temperature at each interior point of the star. For small val- 
ues of (3, the temperature profile is similar to that of the non-causal theory; but 
as (3 is increased, i.e. as relaxational effects grow, it is clear from Figure 1 that 
the temperature profile can deviate substantially from that of the non-causal the- 
ory. Similar results were obtained in the shear-free models studied by Herrera and 
Santos 1 and Govender et al 10 ' 11 . Also, from the plots in Figure 2 the relaxation 
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time for the shear stress exhibits substantially different behaviour when the fluid is 
close to hydrostatic equilibrium as opposed to late-time collapse. In particular, we 
find that 

(Tl)early „ 10Qj (65) 



(n) 



late 



emphasizing the importance of relaxational effects during the different stages of 
collapse. We further note that while the relaxation time for the heat flux is taken 
to be constant, the relaxation time for the shear stress increases as the collapse 
proceeds. 
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Making use of (1) and (52) the proper radius can be written as 




(66) 



Numerical integration of (66) (with Ci(r) = —1) shows that the proper radius is a 
decreasing function of time. This is expected as the star is contracting and losing 
mass. From Figure 3 we note that the proper radius decreases at a faster rate during 
the latter stages of collapse. 

7. Discussion 

We have successfully provided an analytical model of a radiating star undergoing 
gravitational collapse with non-vanishing shear. This model has a Fricdmann-likc 
limit when the heat flux vanishes. We further showed that the causal temperature 
(representing the stellar fluid out of hydrostatic equilibrium) is higher than the 
noncausal temperature at all points of the star. Further analysis revealed that the 
relaxation time for the shear stress (taken to be proportional to the mean collision 
time) increases radially outwards, towards the surface of the star. This is expected, 
as the outer layers of the fluid are cooler than the central regions. Of particular 
significance is the result that the relaxation time for the heat flux (in our case taken 
to be constant) differs from the relaxation time for the shear stress. This is contrary 
to earlier treatments where it was assumed that (r r )^ eat w (T r ) s i iear 5,16 . 

As presented above, the solution (53) together with (54)-(57) admits singulari- 
ties at t — 0,r = and 1 ± Ci(r) exp(3t 1 / 3 /r) = 0. The first singularity is avoided 
by noting that the life of the star is usually taken to start at t = — oo and end at 
t = 0. The model we present here is really suitable for the early life of a star. In 
fact, one can show that a black hole arises as t — ► 0. This consideration also takes 
care of the third case. For early times (large negative values of t in which we only 
take real roots) exp(3t 1 / 3 /r) << 1 and so the expression does not evaluate to zero 
regardless of the sign taken. 

The singularity at r — can be avoided if the solution presented is viewed as 
an "envelope" with a core represented by a metric of the form 



where A and B represent a known static solution (See Govender et aP 7 for a 
further discussion of this approach.). It is a simple matter 28 to match (67) to (53) 
at some inner boundary r . Since (67) is not singular at r = 0, this singularity is 
avoided. 

To make a more realistic comparison of the relaxation times, one requires an 
analytic solution of the causal temperature equation for non-constant relaxation 
times. Future work in this direction will also require the comparison of the various 
relaxation times using the truncated and full transport equations. 



ds 2 = -A (r) 2 dt 2 + B a (r) 2 (f(t) 2 dr 2 + r 2 g(t) 2 dn 2 ), 



(67) 
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